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ABSTRACT. An analysis of the spider assemblage structure and the presence of indicator species in 
the Flemish coastal dunes are presented. The analysis is based on data from more than 170 year-round 
pitfall sampling campaigns from the 1970s onwards. We were able to find indicator species for all identified 
habitats. The assemblages are determined by variation in vegetation structure (succession), atmospheric 
and soil humidity and the occurrence of both natural of anthropogenic disturbance. In the fragmented 
habitats (grasslands and grey dunes), a clear relationship was found between the mean habitat size and 
the stability of the assemblage composition. In moss dominated dunes and short grasslands total species 
numbers do not increase with patch size. Due to microhabitat variation and the possibility of attaining 
viable population sizes the total number of typical species is, however, higher in larger patches. In small 
patches, edge effects are more important and the number of observed species is enlarged by the intrusion 
of species from nearby habitats. 
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Coastal dunes in Europe have been the sub¬ 
ject of several spider community studies 
which reveal that the species composition is 
characterized by the presence of many rare 
and threatened species. Especially in northern 
(Almquist 1973) and western Europe (Duffey 
1968; Bell et al. 1998), such investigations 
have been carried out. Carabid beetles (De¬ 
sender et al. 1992; Desender 1996) and doli- 
chopodid flies (Pollet & Grootaert, 1996) were 
also studied in the Belgian coastal dunes. 
Such studies are necessary for the assessment 
of the conservation value of these habitats. 
General assemblage descriptions together with 
more detailed knowledge of landscape-level 
ecological relationships such as multi-habitat 
use (Bonte et al. 2000a), colonization abilities 
(Bonte et al. 1998), population genetics (De¬ 
sender et al. 1998) and population dynamics 
(Desender 1996; Baert & Desender 1993) 
should be taken into account when developing 
a nature conservation policy. 

Since the beginning of the 20 th century the 
total area of coastal dunes that have not been 


built upon in Belgium diminished from ap¬ 
proximately 6000 ha to less than 3800 ha 
(Vermeersch 1986). The remaining dune areas 
are characterized by an overall increase of 
competitive plant species like Sea Buckthorn 
Hippophae rhamnoides, Burnet Rose Rosa 
pimpinellifolia and Wood Small Reed Cala- 
magrostis epigejos, due to the retreat of local 
dune farmers after World War II and a de¬ 
crease in the rabbit population due to myxo¬ 
matosis and other diseases. This shrub and 
grass encroachment is possibly triggered by 
atmospheric N-deposition and enhanced by 
positive feedbacks in the nitrogen cycle (Veer 
1997) and by the increase of nitrogen-fixing 
Sea Buckthorn. The soil nitrogen and mineral 
content will influence the trophic status of the 
vegetation, which is strongly related to the 
amount of organic components in the upper 
soil layer (Krabbenborg et al. 1983). Current 
habitat management is directed to the conser¬ 
vation and restoration of wet, herbaceous 
grasslands in dune valleys and stable meso- 
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Table 1.—Characterization and total number of pitfall-data from the sampled dune habitats. 


Number of 


Type Indicative plant species pitfall-traps 

Dune woodland Trees: A Inns glutinosa, Acer pseudoplatanus 3 

High, woody shrubs Shrubs: Crataegus monogyna, Hippophae rham- 13 

noides 

Thick humid Calamagrostis grassland Grass: Calamagrostis epigejos 20 

Vital humid Sea-buckthorn-Liguster Shrubs: Hippophae rhamnoides, Salix repens, Li- 10 

shrubs gustrum vulgare 

Wet eutrophic open dune valleys Sedges and grasses: Juncus subnodulosus. Car ex 6 

riparia, Iris pseudacorus 

Thick dry Arrhenaterium grassland Grasses: Arrhenatherium elatius, Avenula pubes- 13 

cens 

Dry Sea buckthorn shrub (in grassland Shrub: Hippophae rhamnoides 6 

mosaics) 

Dwarf shrubs Dwarf-shrub: Rosa pimpinellifolia 18 

Wet mesotrophic open dune valleys Sedges and grasses: Juncus subnodulosus, Car ex 12 

trinervis, C. flacca 

Short grazed mesophytic grasslands Grasses and herbs: Luzula campestris, Galium 15 

verum, Avenula pubescens, Koeleria albescens 

Wet oligotrophic dune valleys Grasses and sedges: Juncus articulatus, Carex tri- 10 

nervis 

Marram dunes Grass: Ammophila arenaria 9 

Moss dominated dry oligotrophic dunes Mosses, annual herbs and grasses: Tortula ruralis 29 

(Grey dunes) ruraliformis, Aira praecox, Erodium circutar- 

ium, Corynephorus canescens 

Bare sand dunes Grass: Festuca rubra arenaria 9 

Anthropogenic distributed sand dunes Herb: Cirsium arvense 5 


phytic grasslands through large-scale removal 
of shrub, followed by horse and cattle grazing. 

In order to develop tools for future moni¬ 
toring we: 1. investigate which parameters in¬ 
fluence the variation in species composition 
and 2. determine indicator spider species for 
all of the dune habitats occurring along our 
coast; and since a crucial question in habitat 
restoration is the effect of patch area on the 
presence of typical species (stenotopic spe¬ 
cies) we also 3. investigate the species-area 
relationships for two highly fragmented hab¬ 
itat types. 

METHODS 

Data collection.—The total assemblage 
analysis is based on data from 178 pitfalls, 
which were operative during an entire year- 
cycle in all kinds of dune habitats of the Bel¬ 
gian coastal dunes from the 1970s onwards 
(Huble 1975; Huble 1976; Van Biervliet 1978; 
Huble & Maelfait 1981; Baert & Desender 
1993; Maelfait 1993; Bonte & Hendrickx 
1997; Bonte et al. 1999; Baert et al. unpub. 
data). In each sampling station three to five 


traps were placed with a distance of 5—10 me¬ 
ter between each pitfall (the traps are glass 
jam jars with a diameter of 9.5 cm, filled with 
a 10% formaline solution). In total more than 
65,000 adult spiders were identified, resulting 
in data on the occurrence of 214 species. 
Voucher specimens are deposited at the Royal 
Belgian Institute of Natural Sciences in Brus¬ 
sels. Of these, 159 species were represented 
by more than five individuals caught and can 
thus be considered resident species and not 
rare vagrants (cf. Maelfait & Baert 1988). The 
sampled vegetation types, the dominant plant 
species and the number of pitfall traps used 
are listed in Table 1. 

Community structure and Indicator spe¬ 
cies.—The community-structure is indirectly 
determined via Detrended Correspondence 
Analysis (Hill 1979a) with the data from the 
separate pitfalls. Only the more abundant spe¬ 
cies were taken into account for the ordination 
analysis. This methodology reveals a multi¬ 
dimensional ordering of the samples (here 
traps) based on their species composition sim- 
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Figure 1.—Dichotomous TWINSPAN clustering of the spider composition based on the species pres¬ 
ence/absence data. (For abbreviations of the assemblages see Table 2). 


ilarity in which traps with a similar assem¬ 
blage are ordered closely together, while those 
with a completely different species composi¬ 
tion are positioned far apart. Since habitat 
characteristics were never studied in a stan¬ 
dardized way during the several sampling 
campaigns, only the habitat type and the 
linked biotic and abiotic variables were indi¬ 
rectly used for the analysis of the parameters 
structuring the spider assemblage. Data on the 
stage of vegetation succession, depth of the 
groundwater level and different kinds of dis¬ 
turbance were taken from Provoost & Hoff¬ 
mann (1996), Different kinds of disturbance 
result from natural eolic dynamics (sand over¬ 
flowing) and anthropogenic factors due to dif¬ 
ferent intensive nature management tech¬ 
niques like mowing of the vegetation (once or 
twice a year) and grazing (year-round) for the 
conservation of oligo- and mesotrophic pas¬ 
tures. 

We used a TWINSPAN-clustering (Hill 
1979b) for the determination of the different 
levels of assemblage similarity, based on the 
species composition from each pitfall trap. In 
this way a hierarchical ordering is obtained in 
which samples are dichotomously separated at 
different levels. The separation at the first lev¬ 
el splits the samples in two different and large 
groups with common species. The samples 
from both groups are then again separated 
based on differences in their species compo¬ 
sition. The total dataset is in this way dichot¬ 
omously clustered at five levels in which the 
highest levels represent the most detailed sam¬ 
ple separation. 

Indicator species for all groups at the dif¬ 


ferent levels from the TWINSPAN-classifica- 
tion were determined with the IndVal-method 
(Dufrene & Legendre 1997). With this meth¬ 
odology, an indicator value is calculated for a 
species in each cluster group. The indicator 
value is calculated as: IndValy = A,-,- * B,-,- * 
100, where A,-,- = N individuals,-, / N indivi¬ 
duals; and By ~ N traps,-,- / N sites,-. In this 
formula Ay is a measure of group specificity, 
where N individuals,^ is the mean number of 
individuals of species i across traps of group 
j and N individuals; is the sum of the mean 
numbers of individuals of species i over all 
groups at that level. By is a measure of fidelity, 
where N traps,-, is the number of traps in group 
j where species i is present, while N sites,- is 
the total number of traps in that group. A,-,- is 
maximum when species i is only present in 
group j, while B„ is maximum when species i 
is present in all traps of the group j. A random 
reallocation procedure of traps among groups 
is used to test the significance of IndVal (500 
permutations). This index (IndVal) is thus 
maximal when all individuals of a species are 
found in a single group of traps and when the 
species occurs in all traps of the group. As a 
consequence the maximal indicator value can 
be interpreted as a measure for habitat speci¬ 
ficity. 

Because pitfall data record (species specif¬ 
ic) activities instead of absolute densities, we 
only analyzed our data by presence/absence in 
the ordination, clustering and IndVal calcula¬ 
tion. In this way, bias to different climatic 
conditions between years are eliminated. 

Assemblage stability and species-area re¬ 
lationships. —The mean Euclidean distance 
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between the different axes-scores for pitfalls 
from the same assemblage, as derived from 
the multi-dimensional DCA-ordination, was 
used as a measure for the assemblage insta¬ 
bility. Low Euclidean distances characterize 
traps that sampled analogue species compo¬ 
sition, while traps with a completely different 
species composition have a high Euclidean 
distance (are ordered distantly). This distance 
measure is thus an indication for the species 
composition similarity between traps from the 
same assemblage. High similarities result thus 
in low distances and indicate stable species 
assemblages. 

Due to shrub and grass encroachment, me- 
sophytic short grasslands and humid valley 
habitats became fragmented and diminished in 
area. A decrease in the area can influence the 
assemblage structure: especially for small 
habitat patches, the assemblage stability is ex¬ 
pected to decrease because of the extinction 
of species and the presence of species from 
adjacent habitats. Therefore, we related the 
habitat (in)stability to the average area of the 
different fragmented habitat types as derived 
from digitized vegetation maps (Provoost & 
Bonte, unpub. data). Differences in stability 
between assemblage groups were assessed 
with Analysis of Variance and related to hab¬ 
itat areas with Spearman rank correlations. 

For moss dominated dunes and short me- 
sophytic grasslands (for which we have exact 
data on the area of the sampled patch), we also 
determined the relation between the total 
number of species and the number of indicator 
species, trapped in three pitfalls, and the area 
of the sampled grassland patch. The species 
number-area relationship was analyzed by 
Pearson correlations between the total number 
of species and between the number of indi¬ 
cator species and the area of the sampled hab¬ 
itat patch, separately for moss dominated 
dunes and short dune grasslands. 

RESULTS 

Assemblage structure.—A total of 15 spi¬ 
der species assemblages were characterized by 
the TWINSPAN clustering. The first division 
clearly separates the eutrophic vegetations 
from the meso- and oligotrophic, short grazed 
habitats. The eutrophic assemblages are sep¬ 
arated at the lowest level in dune woodlands, 
shrubs, marshland and dense grasslands. The 
mesotrophic and oligotrophic habitats are sep¬ 


arated within the second group (Fig. 1) in 
short, rabbit-grazed grasslands, mown meso¬ 
trophic dune valleys, moss dominated dunes, 
dynamic Marram-grass (Ammophiia arenaria) 
dominated dunes and bare dunes. Significant 
indicator species per cluster group (Monte 
Carlo permutations; 500 runs) and their indi¬ 
cator value are listed in Table 2. 

The ordination-analysis clearly shows the 
assemblage structure along three relevant 
axes. The first axis (eigenvalue 0.689) sepa¬ 
rates the different samples along a vegetation 
structure-gradient, where dune woodlands are 
plotted on the left, bare sandy habitats on the 
right. The second axis (eigenvalue 0.587) sep¬ 
arates the humid from the dry habitats: dune 
valley vegetations (dune slacks) and Marram 
dunes above and moss dominated dunes (grey 
dunes) below (Fig. 2). Interesting is the higher 
position along the second axis of Marram 
dunes near the seaside in comparison with 
those along the inner dune front. This stresses 
the importance of atmospheric humidity in ad¬ 
dition to soil humidity as the second important 
assemblage structuring parameter. The third 
relevant axis (eigenvalue 0.383; Fig. 2) is as¬ 
sociated with natural (wind in Marram dunes, 
inundations in dune slacks) and anthropogenic 
dynamics (especially habitat management: 
mowing & grazing in short grazed pastures 
and dune slack meadows); all these disturbed 
habitats which are indeed ordered at the lower 
part of the ordination axis. 

Stability of assemblages from oligo- and 
mesotrophic habitats and species-area re¬ 
lationships.—The assemblage stability differs 
between the different distinguished spider as¬ 
semblages from oligo- and mesotrophic hab¬ 
itats (Bare sand (BS), mown eutrophic valleys 
(ESM), Moss dominated dunes (GD), Marram 
dunes (MD), Mesophytic dune slacks (MS), 
Dry mestrophic grasslands (RGD) and inun¬ 
dating mesotrophic grasslands (RGI) (one 
way-ANOVA, F 1>6 = 11.403, P < 0.001). The 
stability is significantly different between the 
assemblage groups BS, ESM, GD, MD and 
the assemblages of MS, RGI, RGD, but does 
not differ within the two groups. Correlation 
with average patch size is nearly significant 
(Spearman R =—0.750, P = 0.052) and in¬ 
dicates that assemblages from small habitats 
tend to be more diverse in species composi¬ 
tion. 

The species-area relationship of the total 
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Table 2.—Indicator species (Monte Carlo permutations, P < 0.01) and indicator value (IndVal) at the 
different cluster levels (See Fig. 1), with description of the assemblage habitat characteristics. 


Habitat (abbreviation) and habitat characteristics 

Indicator species 

IndVal 

Eutrophic vegetation (EV) 

A lope cos a pulverulenta (Clerck, 1757) (Ly- 

74.97 

Higher, litter rich, dense 

cosidae) 


Dry or humid 

Bathyphantes parvulus (Westring, 1851) 
(Linyphiidae) 

51.04 


Centromerus prudens (O.P.-Cambridge, 

1873) (Linyphiidae) 

60.06 


Centromerus sylvaticus (Blackwall, 1841) 
Linyphiidae) 

87.08 


Clubiona comta C.L. Koch, 1839 (Clubioni- 
dae) 

23.60 


Clubiona lutescens Westring, 1851 (Clubion- 
idae) 

25.71 


Episinus angulatus (Blackwall, 1836) (Ther- 
idiidae) 

16.85 


Ero furcata (Villers, 1789) (Mimetidae) 

47.03 


Euryopus fiavomaculata (C.L. Koch, 1836) 
(Theridiidae) 

49.49 


Floronia bucculenta (Clerck, 1757) (Linyphi¬ 
idae) 

27.77 


Gonatium rubens (Blackwall, 1833) (Liny¬ 
phiidae) 

69.57 


Linyphia triangularis (Clerck, 1757) (Liny¬ 
phiidae) 

12.94 


Maso sundevalli (Westring, 1851) (Linyphi¬ 
idae) 

43.76 


Meioneta saxatilis (Blackwall, 1844) (Liny¬ 
phiidae) 

46.09 


Neriene clathrata (Sundevall, 1830) (Liny¬ 
phiidae) 

25.97 


Ozyptila simplex (O.P.-Cambridge, 1862) 
(Thomisidae) 

71.09 


Palliduphantes ericaeus (Blackwall, 1853) 
(Linyphiidae) 

19.97 


Palliduphantes pallidus (O.P.-Cambridge, 
1871 (Linyphiidae) 

69.56 


Pirata hygrophillus Thorell, 1872 (Lycosi- 
dae) 

49.16 


Pocadicnemis juncea Locket & Millidge, 
1953 (Linyphiidae) 

59.63 


Robertas lividus (Blackwall, 1836) (Theridi¬ 
idae) 

74.39 


Theridion bimaculatum (Linnaeus, 1758) 
(Theridiidae) 

33.35 


Trochosa terricola Thorell, 1856 (Lycosidae) 

88.01 


Walckenaeria acuminata Blackwall, 1833 
(Linyphiidae) 

47.11 


Walckenaeria atrotibialis (O.P.-Cambridge, 
1878) (Linyphiidae) 

93.17 


Zora spinimana (Sundevall, 1833) (Zoridae) 

66.94 

Meso-oligotrophic vegetation (MOV) 

Arctosa perita (Latreille, 1799) (Lycosidae) 

63.24 

Short, sparse vegetation 

Haplodrassus dalmatensis (L. Koch, 1866) 

73.01 

Rabbit grazed 

(Gnaphosidae) 


Dry or humid 

Meioneta rurestris (C.L. Koch, 1836) (Liny¬ 

34.49 

Sandy patches 

phiidae) 
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Table 2.—Continued. 


Habitat (abbreviation) and habitat characteristics Indicator species IndVal 


Woodland and woody shrubs (WS) 

High vegetation and litter rich 
Presence of trees {Crataegus monogyna) 

Dense grasslands (EG) 

Dense and tall grass layer 

Litter-rich 

Dry or humid 


Thermophilious grasslands (TG) 
Short, no or scarce litter 
Dry or humid 

Dynamics: wind, grazing or mowing 


Bare sand (BS) 

No vegetation 
Dune woodland (DW) 

Dominance of trees (Ainus glutinosa) 

Litter-rich 

Humid 


Woody shrubs (SC) 

Dominance of Sea Buckthorn {Hippophae rham- 
noides) and Calamagrostis epigejos 
Presence of trees {Crataegus monogyna) 


Parapelecopsis nemoralis (O.P.-Cambridge, 53.85 
1884) (Linyphiidae) 

Styloctetor romanus (O.P.-Cambridge, 1872) 37.22 

(Linyphiidae) 

Tegenaria agrestis (Walckenaer, 1802) (Age- 24.72 
lenidae) 

Xysticus sabulosus (Hahn, 1832) (Thomisi- 40.28 
dae) 

Tapinopa longidens (Wider, 1834) (Linyphi- 19.35 
idae) 

Walckenaeria nudipalpis (Westring, 1851) 44.13 

(Linyphiidae) 

Clubiona diversa O.P.-Cambridge, 1862 8.09 

(Clubionidae) 

Cnephalocotes obscurus (Blackwall, 1834) 35.07 

(Linyphiidae) 

Enoplognatha thoracica (Hahn, 1833) (Ther- 30.89 
idiidae) 

Ero cambridgei Kulczynski, 1911 (Mined- 11.53 

dae) 

Pachygnatha degeeri Sundevall, 1830 (Te- 71.92 
tragnathidae) 

Pisaura mirabilis (Clerck, 1757) (Pisauridae) 18.92 

Walckenanaeria antica (Wider, 1834) (Liny- 41.61 
phiidae) 

Agroeca lusatica (L. Koch, 1875) (Liocrani- 37.57 
dae) 

Alopecosa barbipes (Sundevall, 1833) (Ly- 34.59 
cosidae) 

Bolyphantes luteolus (Blackwall, 1833) (Lin- 16.00 
yphiidae) 

Centromerita concinna (Thorell, 1875) (Lin- 62.25 
yphiidae) 

Walckenaeria monoceros (Wider, 1834) (Lin- 32.82 
yphiidae) 

Xysticus kochi Thorell 1872 (Thomisidae) 64.91 

Erigone longipalpis (Sundevall, 1830) (Lin- 57.14 
yphiidae) 

Ceratinella scabrosa (O.P.-Cambridge, 1863) 43.48 

(Linyphiidae) 

Diplocephalus picinus (Blackwall, 1841) 82.86 

(Linyphiidae) 

Enoplognatha ovata (Clerck, 1757) (Theri- 14.29 
diidae) 

Macrargus rufus (Wider, 1830) (Linyphiidae) 19.05 

Pardos a saltans Topfer-Hofmann, 2000 (Ly- 66.67 
cosidae) 

Tapinocyba insecta (L. Koch, 1869) (Liny- 56.60 
phiidae) 

Tenuiphantes zimmermanni (Betkau, 1890) 85.19 

(Linyphiidae) 

Monocephalus fuscipes (Blackwall, 1836) 88.92 

(Linyphiidae) 

Saaristoa abnormis (Blackwall, 1841) (Lin- 44.44 
yphiidae) 
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Table 2,—Continued. 


Habitat (abbreviation) and habitat characteristics Indicator species IndVal 


Eutrophic wet dune valleys (ES) 
Humid, Winter inundations 

High, dense vegetation 
Dominance of Car ex riparia 


Mesotrophic grasslands (SG) 

Dry or humid (winter indundations) 
Marram and moss dominated dunes (MG) 
Sandy, scarce vegetation 
Mainly mosses and lichens 
Ammophila arenaria- tussocks 


Antropogenic disturbed sandy soils (SA) 

Bare sand, human activities 
Dense shrubs (DS) 

Dominance of Hippophae rliamnoides and Ligus- 
trum vulgare 


Degradating Shrub (CS) 

Shrub with open patches, colonized by Calama- 

grostis epigejos 


Centromerita hicolor (Blackwall, 1833) (Lin- 42.17 
yphiidae) 

Ceratinella brevipes (Westring, 1851) (Liny- 29.78 
phiidae) 

Clubiona reclusa O.P.-Cambridge, 1863 34.38 

(Clubionidae) 

Dicymbium nigrum (Blackwall, 1834) (Liny- 63.62 
phiidae) 

Gnathonarium den tat um (Wider, 1834) (Lin- 15.63 
yphiidae) 

Pachygnatha clercki Sundevall, 1823 (Te- 50.86 

tragnathidae) 

Pardos a palustris (Linnaeus, 1758) (Lycosi- 83.92 
dae) 

Pardosa pullata (Clerck, 1757) (Lycosidae) 71.79 

Pirata latitans (Blackwall, 1849) (Lycosidae) 85.49 

Pirata piraticus (Clerck, 1757) (Lycosidae) 42.19 

Tiso vagans (Blackwall, 1834) (Linyphiidae) 71.12 

Troxochrus cirrifrons (O.P.-Cambridge, 30.95 

1871) (Linyphiidae) 

Troxochrus scabrosa (Westring, 1851) (Lin- 34.57 
yphiidae) 

Pardosa monticola (Clerck, 1757) (Lycosi- 53.43 
dae) 

Agroeca cnprea Menge, 1873 (Liocranidae) 66.98 

Drassodes cupreus (Blackwall, 1834) (Gna- 50.05 
phosidae) 

Dysdera croc at a C.L. Koch, 1838 (Dysderi- 17.81 
dae) 

Metopobactrus prominulus (O.R-Cambridge, 37.83 

1872) (Linyphiidae) 

Poeciloneta variegata (Blackwall, 1841) 15.79 

(Linyphiidae) 

Sitticus saltator (O.P.-Cambridge, 1868) (Sal- 32.50 
ticidae) 

Thanatus striatus C.L. Koch, 1845 (Thomis- 39.97 
idae) 

Pardosa proxima (C.L. Koch, 1847) 57.14 

(Lycosidae) 

Agyneta subtilis (O.P.-Cambridge, 1847) 29.70 

(Linyphiidae) 

Gongylidium mfipes (Linnaeus, 1758) (Lin- 61.78 
yphiidae) 

Microneta varia (Blackwall, 1841) (Linyphi- 24.48 
idae) 

Ozyptila praticola (C.L. Koch, 1837) (Thom- 31.43 
isidae) 

Pholcomma gibbum (Westring, 1851) (Ther- 11.54 
idiidae) 

Walckenaeria cucculata (C.L. Koch, 1836) 43.50 

(Linyphiidae) 

Agyneta decora (O.P.-Cambridge, 1871) (Lin- 10.77 
yphiidae) 
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Table 2.—Continued. 


Habitat (abbreviation) and habitat characteristics 
Humid 


Wet rough litter rich vegetation (ESR) 

Rough, eutrophic vegetation 

Inundations, no management 

Dry dense grasslands-shrub mosaics (AG) 

Mosaics of low shrubs and Avem</c/-grassland 

Dry, no management 


High dwarf shrubs (RP) 

Dominance of Rosa pimpinellifolia and Arrhena- 
terium elatius 
High grass layer 
Presence of litter 


Mesotrophic dune valleys (MS) 
Dominance of Junciis subnodulosiis 
Yearly mowed 
Winter inundations 


Short Mesoptrophic grasslands (RG) 

Wet (inundating) or dry, rabbit grazed 
Marram dunes (MD) 

Dominance of Marram grass ( Ammophila 
arenaria ) 

Strong wind dynamics, close to the sea 
Scarce vegetation 


Indicator species IndVal 

Ceratinella brevis (Wider, 1834) (Linyphi- 22.00 
idae) 

Kaestneria pullata (O.P.-Cambridge, 1863) 41.95 

(Linyphiidae) 

Clubiona phragmites C.L. Koch, 1843 (Clu- 35.71 
bionidae) 

Xysticus ulmi (Hahn, 1831) (Thomisidae) 16.67 

Agyneta conigera (O.P.-Cambridge, 1863) 11.48 

(Linyphiidae) 

Hahnia nava (Blackwall, 1841) (Hahniidae) 31.79 


Metellina mengei (Blackwall, 1870) (Tetrag- 14.29 
nathidae 

Maso gallicus Simon, 1894 (Linyphiidae) 20.21 
Philodromus cespitum (Walckenaer, 1802) 17.54 

Alopecosa cuneata (Clerck, 1757) (Lycosi- 28.18 


dae) 

Heliopliamis flavipes (Hahn, 1832) (Saltici- 26.67 
dae) 

Xysticus erraticus (Blackwall, 1834) (Thom- 36.42 
isidae) 

Trachyzelotes pedestris (C. L. Koch, 1837) 27.78 

(Gnaphosidae) 

Arctosa leopardus (Sundevall, 1833) (Lycos- 33.79 
idae) 

Clubiona trivialis C.L. Koch, 1843 (Clubion- 36.11 
idae) 

Collinsia innerans (O.P.-Cambridge, 1855) 20.00 

(Linyphiidae) 

Erigone arctica (White, 1852) (Linyphiidae) 65.51 

Erigone promiscua (O.P.-Cambridge, 1873) 74.74 

(Linyphiidae) 

Prinerigone vagans (Audouin, 1826) (Liny- 56.45 
phi idae) 

Gongylidiellum vivum (O.P.-Cambridge, 39.46 

1975) (Linyphiidae) 

Oedothorax apicatus (Blackwall, 1850) (Lin- 41.59 
yphiidae) 

Oedothorax fuscus (Blackwall, 1834) (Liny- 74.22 
phiidae) 

Oedothorax retusus (Westring, 1851) (Liny- 58.23 
phiidae) 

Thyphochrestus digitatus (O.P.-Cambridge, 42.21 
1872) (Linyphiidae) 

Clubiona frisia Wunderlich & Schuett, 1995 75.41 
(Clubionidae) 

Clubiona subtilis L. Koch 1867 (Clubioni- 31.07 
dae) 

Micarici pulicaria (Sundevall, 1831) (Gna- 55.56 
phosidae) 

Porrhomma microphthalmum (O.P.-Cam- 14.81 
bridge, 1871) (Linyphiidae) 

Tibellus maritirnus (Menge, 1875) (Thomisi- 18.52 
dae) 

Trochosa ruricola (De Geer, 1778) (Lycosi- 
dae) 


28.70 
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Table 2.—Continued. 


Habitat (abbreviation) and habitat characteristics Indicator species IndVal 


Moss dominated dunes & Marram dunes near the 
inner dune front (GD) 

Dominance of lichens and mosses 
Scarce Marram grass vegetation 
Near inner dune front 


Dry mesotropic grasslands (RGD) 

Rabbit grazed, short grass layer 
Dominance of Luzula campestris 

Inundating mesotrophic grasslands (RGI) 
Inundating, short Career-vegetation 
Presence of Creeping willow {Salix repens) 


Alopecosa fabrilis (Clerck, 1757) (Lycosidae) 20.85 

Micaria dives (Lucas, 1846) (Gnaphosidae) 18.97 

Walckenaeria stylifrons (O.P.-Cambridge, 

1875) (Linyphiidae) 30.38 

Zelotes longipes (L. Koch, 1866) (Gnaphos- 66.70 
idae) 

Pelecopsis parallella (Wider, 1834) (Liny- 35.81 
phiidae) 

Trichopterna cito (O.P.-Cambridge, 1872) 45.98 

(Linyphiidae) 

Cheiracanthium vi re see ns (Sundevall, 1833) 51.43 

(Clubionidae) 

Lepthothrix hardyi (Blackwall, 1850) (Liny- 26.98 
phiidae) 

Xeroiycosa miniata (C.L. Koch, 1834) (Ly- 64.55 
cosidae) 


number of species and the total number of 
specific (indicator) species as a function of the 
area of moss dominated dune and short dune 
grassland patches is illustrated in Fig. 3. The 
relationship between patch size and total num¬ 
ber of species is not significant for either the 
moss dominated or the short dune grasslands 
(Fig. 4: Pearson correlation, r < 0.20; P > 
0.05). The number of resident indicator spe¬ 
cies however is higher in large patches in both 
vegetation types (Fig. 4: Pearson correlation 
for moss dunes: r = 0.87, P < 0.01 and for 
dry mesoptrophic dune grasslands: r = 0.93, 
P < 0,01). 

DISCUSSION 

Our results indicate that almost ail dune 
system habitat types are characterized by the 
presence of indicator species, dependent on 
the cluster level. Desender (1996) showed that 
typical dune carabid species have a strong 
year-to-year fluctuation in population size, but 
were never completely lacking from the sam¬ 
ples. This variation could be explained by var¬ 
iation in climatological variables. We there¬ 
fore used only absence/presence data, so true 
indicators that are always present (indepen¬ 
dent of their yearly abundance) are unambig¬ 
uously identified. Besides year-to-year fluctu¬ 
ations, species assemblages can vary as a 
function of habitat conditions and landscape 
structure. Our analysis is based on an exten¬ 
sive data set from habitats of different size and 


from different landscape configurations, so the 
determined indicator species can be used as 
bio-indicators for future monitoring of the 
management of both open (dominance of 
grasslands, Marram dunes) and closed (shrub 
dominated) dune landscapes. 

That fact that the variation in species as¬ 
semblages can be explained by variation in the 
vegetation structure or succession stage is not 
surprising and has been documented several 
times in other studies of invertebrate assem¬ 
blages (spiders: Duffey 1968; Almquist 1973; 
carabid beetles: Desender et al. 1992; Empidid 
and dolichopdid flies: Pollet & Grootaert 
1996). Our study indicates the importance of 
atmospheric and soil humidity as a second im¬ 
portant assemblage structuring component 
since wet vegetation types are clearly sepa¬ 
rated from dry ones, and Marram dunes near 
the seaside were separated from those of drier 
more inland dunes. The importance of atmo¬ 
spheric humidity is demonstrated by the pres¬ 
ence of species from dune valleys (Clubiona 
/m/a and C, subtilis) in Marram dunes near 
the seaside, while they are completely absent 
from the same habitat near the inner dune 
front (1-2 km from the seaside), where at¬ 
mospheric humidity is significantly lower 
(Provoost & Hoffmann 1996). The same phe¬ 
nomenon (defined as a double ecological oc¬ 
currence) has also been documented by Duf¬ 
fey (1968) in British coastal dimes. A third 
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Figure 2.—DCA ordination of the pitfall data based on the species absence/presence data with indication 
of the habitat type. Above: ordination along the first and second axes; below: ordination along the first 
and third axes. 


important abiotic source of variation is de¬ 
fined here as habitat disturbance. The third 
axis separates stable habitats without distur¬ 
bance (woodland, shrubs, dwarf shrubs, rough 
permanent grasslands) from habitats with high 
natural (inundations: wet open dune slacks; 


eolic dynamics: Marram dunes, bare dunes) or 
anthropogenic disturbance (grazing and mow¬ 
ing management). These habitats are charac¬ 
terized by ruderal species like Erigone atra, 
E. dentipalpis, E. arctica, Oedothorax fuscus, 
O, retusus, O. apicatus and Bathyphantes 
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Figure 3.—Assemblage instability of fragmented grassland and dune valley habitats (mean Euclidean 
distances of DC A scores and 95% confidence intervals; n = number of distance measurements) in function 
to the average habitat patch size (means and 95% confidence intervals). For abbreviations see Table 2. 



gracilis which are all short-living spiders with 
a rapid juvenile development and a well de¬ 
veloped ballooning dispersal capacity. 

This data analysis is based only on the pres¬ 
ence of adult spiders. Earlier studies have in¬ 
dicated that species typical for open habitats 
like short grasslands and moss dominated dry 
dunes (grey dunes) need proximate patches of 
dense and litter-rich vegetation for their ju¬ 
venile development and/or retreat during un¬ 
favorable periods in their mature life-stage 
(Bonte et al. 2000 a & b). Thus, habitat var¬ 
iation can strongly alter the presence of spider 
species bound to these dense vegetation patch¬ 
es for their juvenile development. Although 
not documented for spider populations, a min¬ 
imal patch area can determine the presence of 
viable population size (Hanski 1999). In both 
cases, an increasing patch area should affect 
the spider assemblage directly because patch 
area influences the population size or indi¬ 
rectly because an increasing patch size en¬ 
hances internal microhabitat variation. Our re¬ 
sults on the assemblage level show that the 


stability in species composition of spider as¬ 
semblages in patchy habitats depends on the 
mean patch area, indicating that differences in 
the spider assemblage vary more in small hab¬ 
itat patches than in larger ones. Edge effects 
in small habitats can alter the spider assem¬ 
blage dramatically, because of the intrusion of 
species typical for other habitat types in the 
patch, due to the high circumference-surface 
ratio. This is certainly true for moss dominat¬ 
ed dunes and short grasslands: total species 
numbers do not differ as a function of the 
patch size while the number of indicator spe¬ 
cies significantly increases with an increasing 
patch size. An explanation of this species-area 
relationship cannot be given without further 
research on both internal micohabitat hetero¬ 
geneity and minimal population sizes. Varia¬ 
tion in soil conditions can explain the aggre¬ 
gation of soil-dwelling arthropods like 
springtails Collembola (Bonte & Mertens, un¬ 
pub. data). Since these are the main prey for 
typical juvenile wolf spiders and adult dwarf 
spiders, a larger patch size can alter the total 
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Ln(area) (m 2 ) 

Figure 4.—Species-area relationship for all species and indicator species caught on moss dominated 
dry dunes and short grasslandsfrom different sizes (Spearman corelation, *: P < 0.01). 


number of indicator species indirectly by the 
presence of higher internal microhabitat vari¬ 
ation. For the study of minimal patch size and 
related population size, more detailed studies 
are needed on meta-population dynamics, 
based on the survey of a higher number of 
habitat patches. 
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